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Модификация метода Литтла для решения 
кольцевой задачи о сельском почтальоне 


В статье приведена формулировка и дан анализ обобщения гамильтоновой задачи о сельском почтальоне. 
Предлагается точный метод её решения, развивающий классический алгоритм Литтла. 


Введение 


Далеко не для каждой задачи оптимизации на транспортной сети разработаны 
алгоритмы поиска решений, пригодные в реальных ситуациях. Одной из таких задач 
является задача о сельском почтальоне, ограниченная версия которой рассматри- 
вается в данной статье. Эту версию можно рассматривать как гамильтонову задачу 
коммивояжера (ГЗК) с условием, сужающим пространство её решений. 


Постановка задачи. Задан связный взвешенный граф Н = (7,0) с множеством 
вершин Г, |ИИ=и и множеством рёбер И. Каждому ребру {1,ЛеОИ приписан вес 
Ч; Е 70, 1%], 5] =Ьп, 75 — множество неотрицательных целых чисел. Граф Н 


полностью определяется симметричной матрицей стоимостей [4 | ‚ где а; = . 


И 


ли 1 ДЛЕЦи а. =0, иначе 17, т] =Ьп, а. =0, 1=Ь п. мн ТВ 
ес Чи 4, =® аче 1= 1 4; = 1 На множестве И 


задано непустое подмножество рёбер А<И. 

Требуется найти в графе Н простой цикл, включающий каждое ребро А и 
имеющий минимальную сумму весов всех рёбер. 

В рассматриваемой задаче в отличие от задачи о сельском почтальоне (3СП) 
искомый цикл должен быть простым. Поставленную задачу назовём кольцевой зада- 
чей о сельском почтальоне (КЗСП). 


Пусть У(К) — простой цикл графа Н, содержащий все рёбра заданного мно- 


жества А. Тогда КЗСП состоит в нахождении простого цикла у’ (К), доставляющего 


С(›(в))= У 4. я 


|, Ве) 


минимум функционалу 


Основная идея метода 


Ясно, что КЗСП М№Р-полна, и только её очень ограниченные частные случаи 
эффективно разрешимы. 

Из М№Р-полноты КЗСП следует, что она не поддается эффективным точным ме- 
тодам решения, а из того, что множество простых циклов графа, включающих все 
рёбра А, может оказаться пустым, вытекает неприменимость эффективных прибли- 
женных алгоритмов минимизации (1). Таким образом, единственной приемлемой 
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схемой поиска решения КЗСП является схема сокращения перебора, организованная 
по методу ветвей и границ, который остаётся универсальным и мощным инструмен- 
том комбинаторной оптимизации. 


В представленном здесь методе поиска цикла у’(К) алгоритму ветвей и гра- 


ниц предшествует стадия эффективной проверки известных достаточных условий 
неразрешимости КЗСП. К ней относится процедура вершинно-рёберного преобразо- 
вания (ВРП) исходного графа Н в граф со структурными характеристиками, вызы- 
вающими выполнение переборного алгоритма [1]. 

Переборный алгоритм, выполняемый на второй стадии решения КЗСП, пред- 
ставляет собой модификацию классического метода Литтла [2], адаптированного для 


нахождения цикла у“ (В) связного подграфа Н = (У, И) исходного подграфа, в кото- 


ром: а) степень каждой вершины ТЕТ больше 2, 6) множество А образует паросоче- 
тание, в) нет точек сочленения. 

В рассматриваемой модификации, как и в алгоритме Литтла для решения сим- 
метричной задачи коммивояжера (ЗК), дерево ветвлений развивается из представ- 


ления матрицы стоимостей | 4. | ориентированным графом С = (У, Е), полученным 
й 


в данном случае из графа Н заменой каждого ребра {1, }е И парой дуг (1, ЛЕЁ и 
(7, ЛЕЕ. Решением КЗСП является простой контур минимальной стоимости. Следо- 
вательно, модификация алгоритма применима для нахождения решения КЗСП, сфор- 
мулированной в терминах ориентированных графов. 

В методе решения поставленной задачи корневой вершине Хо дерева перебора 


ставится в соответствие матрица стоимостей О = Га, исходного графа Н, матрица 
ь п 


длин кратчайших путей М и матрица маршрутизации И’. После их преобразования 
находятся вершины, порождённые на шаге ветвления. Для определения матриц М и 
Й’ на каждом шаге ветвления применяется модификация алгоритма Флойда-Уор- 
шелла. В совокупности, матрицы О, М и И’ позволяют выполнять в условиях постав- 
ленной задачи все действия, включенные в классический метод Литтла. В приве- 
денной версии метода выбор очередной вершины ветвления выполняется по схеме, 
изображенной на рис. 1. 

Из матриц М и ПИ’, формируемых модифицированным алгоритмом Флойда- 
Уоршелла для текущих параметров матрицы О, определяется дуга (р, 4) или путь из 
вершины р в вершину 4, инициирующие ветвление. 

Условия рассматриваемой задачи требуют формирования множества О запре- 
щенных дуг, т.е. дуг, которые приводят к появлению подконтуров в искомом реше- 
нии у“ (К) : 

В множество О включаются дуги в зависимости от способа разбиения прост- 
ранства допустимых решений на подмножества. 

Если ветвление инициирует дуга (А, /), соответствующая ребру {^, ПЕК, то 
множество допустимых решений задачи разбивается на два подмножества. Первое 
подмножество состоит их всех контуров, включающих дугу ((, [), а второе — из всех 
контуров, содержащих дугу ([, А). 

Пусть ветвление вызывает дуга (А, [), {К, 1 е К, или путь из вершины К в вер- 
шину /. При формировании множества О он рассматривается как дуга ((, /). В этом слу- 
чае множество допустимых решений разбивается на два подмножества. Первое множест- 


во включает все контуры, содержащие ((, /), а второе — не содержащие ((, /), т.е. (А,!). 
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Р- матрица стоимостей 
Модифицированный алгоритм Флойда-Уоршелла 


М - матрица длин Й/ — матрица 
кратчайших путей маршрутизации 


Выбор дуги или пути, 


Мьь - рабочая матрица 
инициирующего ветвление 


Процедура 
ветвления 


Рисунок 1 — Схема ветвления в методе решения КЗСП 


При формировании первого подмножества требуется проверка следующих усло- 
вий. Если множество А содержит ребро {х, А}, то в искомый контур вместе с дугой 
(К, [) включается дуга (х, К). Аналогично, если множество К содержит ребро {у, 1}, то 
к дуге (К, /) присоединяется дуга (1, у). Включение дуги (х, К) или (ру) в подмно- 
жество решений означает, что определяющая его матрица О не содержит не только 
строку К и столбец [, но и ту строку и столбец, номера которых являются началом и 
концом присоединенной дуги. В случае {х, А}, {у, ПЕК кдуге (А, [) присоединяются 
дуги (х, К) и (ру), а в соответствующей матрице О, определяющей это множество, 
исключаются строки х, А, [и столбцы (, [, у. 

Сформулируем правила определения множества О: 

— для корня Хо дерева перебора О =@; 

— все элементы множества О при вершине ветвления передаются её порожден- 
ной вершине; 

— если в любое решение при вершине, порожденной на шаге ветвления, вклю- 
чается дуга (А, /), то к множеству О при этой вершине добавляется дуга (1, К); 

— вес каждой дуги в О принимается равным бесконечности. 

Если ветвление вызывает путь из вершины К в вершину /, включающий не ме- 
нее двух дуг, то для разбиения множества допустимых решений на непересекающи- 
еся подмножества используется правило, которое отличается от правила, применяе- 
мого для формирования множества 0. В этом случае результатом разбиения является 
множество вершин И(Ё), представляющее расширение множества активных вершин, 
порожденных на предыдущих шагах ветвления. 

В корневой вершине и вершине ветвления, для которой процесс разбиения ини- 
циирует дуга, /(Г.) = ©. Множество вершин ИК) ИР) определяет рабочую под- 
матрицу матрицы М. Рабочая матрица не содержит строк и столбцов, соответству- 
ющих вершинам множества /\(ИК)‹ ИГ.)). В результате приведения рабочей матрицы 
находится рабочая приведенная матрица Мег. 
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Пусть ветвление инициирует путь, содержащий несколько дуг. Тогда началь- 
ная и конечная вершины пути принадлежат множеству ИК) ИГ), а его промежу- 
точные вершины — множеству Г/\(ИЮ) КЕ.)). Действительно, начальная вершина 
соответствует строке рабочей матрицы, конечная вершина — столбцу, а рабочая 
матрица построена на вершинах множества ИА)» ИЕ). Промежуточная вершина 


У,,1=2, &-—1 пути (› ре ) требует включения соответствующей ей 


Я › Ур а УЖ ь 


строки и столбиа в рабочую матрицу. В дереве перебора вершине у, соответствует 


множество решений, содержащих эту вершину. Каждое такое множество можно рас- 
сматривать как разбиение на два подмножества, одно из которых содержит дугу 


(> У т ‚ а другое не включает её. В дереве перебора путь (> ) порож- 


дает А висячих вершин ХльХ 2» т а. . 


ти Х, расширение совпадает с расширением ИГ) для Х я за при вер- 


Л? У > 2 УЖ 
исходящих из вершины Х ГЕ При 


вершинах Х 


шине Х,,,/[=2, &-1 — оно равно ИВ, |. 


> 
Матрицы О при висячих вершинах формируются из матрицы О при вершине 
ветвления Х я следующим образом: 


1. В матрице ДР для вершины Х„;,{=1,К-1 элемент ( Е На) принимается рав- 


НЫМ 00. 


2. Чтобы получить матрицу О для вершины Х д › ИЗ матрицы для Х я исключа- 


ются те строки и столбцы, по которым строится путь (> НУ ) | 

3. Матрица Р при вершине Х, 1=2, К-1 содержит строку, соответствующую 
вершине у, ‚ ине содержит строк и столбцов, определяющих путь (> Ир» ++ Уй ) ь 

Множество ИР) формируется согласно следующему правилу: 

— в корневой вершине ИГ) =©; 

— все элементы множества ИГ) передаются вершине, порожденной вершиной 
ветвления; 


— если ветвление вызывает путь (5, ...,р, ..., 9, включающий А промежуточных 
вершин, то соответствующее ему множество решений с расширением ИГ.) разбива- 
ется на последовательность из К + 2 подмножеств ($, ..., Р,..., Г) с расширением 
ИГ) {р} для подмножества Р, р=5, #; в любом решении подмножества Р содер- 
жится часть пути (5, ...,р, ..., В) из вершины 5 в вершину р; каждое решение подмно- 
жества Т включает путь (5, ..., р, ..., ); все решения подмножества 5 не включают 
этого пути. 


Модификация алгоритма Флойда-Уоршелла 


Платформой, на которую опираются полученные правила ветвления, служит 
модификация известного алгоритма Флойда-Уоршелла, определяющего в ориенти- 
рованном графе кратчайшие пути между всеми парами вершин [3]. Модификация 
формирует матрицы М и И’ в каждой точке ветвления в результате выполнения сле- 
дующих действий. 
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Вход. р — матрица стоимостей 4, подграфа С(О) ориентированного графа @ = 
=(У, Е), Е {1,2,..., п}; С(О) = С в корневой вершине дерева перебора; 

К — множество рёбер, образующих паросочетание исходного графа Н; 

О- множество запрещённых дуг в вершине ветвления, 0 =©) в случае С(О) = С; 

ИК) — множество вершин, инцидентных рёбрам паросочетания К; 

ИГ) - расширение множества вершин в точке ветвления; ИГ) = © в корне де- 
рева перебора; 

оу» [со/— множество индексов соответственно строк и столбцов матрицы 2 по- 


рядка | Ги» | = | [со |, п — порядок матрицы ДО в корневой вершине; 
Выход. М — матрица длин т, кратчайших путей (41, К, ..., 1,7), Б/Е ИК), К, 1 ИК), 
ЛЕК; 


Й’ — матрица кратчайших путей у›;, соответствующая М. 
Подготовительный этап: 


1% Гогай ТЕ [лю 90 


я Гог а уе Гсо/ 40 
3°. Бест 
с со, если [= ] или т; = 5, 
4. ее ры. 
У, если! * } ит, < о. 
0 == 
6%. епд 


Этап нахождения элементов матриц Ми И’: 


1. Юга КЕ о! 90 

2 Нке ИК) ИЕ) еп 

3 Гог а] те [ох 90 

4. Гог а| уе /со/ 90 

Э НТ» ] апд 1=К ап@д 7 =К еп Беот 
6 НИ 1, ДЕК еп сопйптие 

7 И (1) ЕО Феп сопйтие 

8 И ту> ти ть {Пеп 


а. 


Бест 
10. ту = ти + ть 
Гр епа 
13. епа 


На каждом шаге ветвления множеству О запрещенных дуг ставится в соответ- 
ствие подграф С(О) ориентированного графа С = (Т, ЕР) и его матрица стоимостей О. 

В строке 2 алгоритма содержится условие, исключающее возможность прохож- 
дения путей через вершины множества ИК) И(Ё), для которого формируется рабо- 
чая матрица М»и. Так как в корне дерева перебора ИГ)=Ф© ‚ то в этом случае Мьи яв- 
ляется квадратной матрицей порядка |И(К)|. 

Строка 6 для каждого ребра {1, /} е К запрещает построение кратчайшего пути 
между вершинами 1 и / и сохраняет веса дуг (1, /) и (7, 1 в матрице О. 

Дуги множества О, запрещенные в процессе ветвления для устранения недо- 
пустимых контуров, должны быть запрещены и при формировании матриц М и И. 
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Как и алгоритм Флойда-Уоршелла, предложенная модификация имеет времен- 
ную сложность, оцениваемую полиномом третьей степени от размерности входной 
матрицы 0. 


Нижние границы стоимости кольцевых маршрутов 


Нижняя оценка стоимости любого решения в данной точке ветвления опреде- 
ляется из подматрицы Маг рабочей матрицы М. В матрицу М, полученную из мат- 
рицы О модифицированным алгоритмом Флойда-Уоршелла, включены все строки с 
индексами множества /„„ и все столбцы с индексами множества /соь, | [мо | = | Го |, 
Трои» [со © {1,2,..., п}. Матрица Миг состоит из тех элементов, которые принадлежат 
строкам /„„„б (ИСК) И) и столбцам и“ (ИКЮ) И )). 

Для нахождения нижней оценки выполняется приведение матрицы Мрг, вклю- 
чающее такие действия: 

1) найти в каждой строке минимальный элемент 2; и вычесть его из каждого эле- 
мента этой строки; 

2) исключить из рассмотрения все строки, соответствующие вершинам множества 
ИГ), в каждом столбце полученной подматрицы найти минимальный элемент Й; и 
вычесть его из каждого элемента этого столбца; 

3) в полученной матрице найти те элементы (1, /) и (1,7, {1 ДЕК, для которых 


Д; = пай ту, т’, й ‚ и положить 


ее р те -АОН | 
т; =т,; Д;,, если т, * 0, т, =ту Д;, если т» *о. 


Пусть Ё, — совокупность всех дуг любого решения из подмножества решений, 
представленного точкой ветвления (. Оценка снизу в точке { определяется по формуле: 


ф = уз 8+ » В+» Д+ 2%. (2) 
оС (В)УТ(Г,)} ЛЕГсо1 (Е) ,ЛеЕЮ (ЕЁ, 
1,1 ЕЕ} 
ЛЕ 
Первые два слагаемых в формуле (2) позволяют оценить стоимость маршрутов 
так, как и в методе Литтла. Третье слагаемое увеличивает нижнюю границу, исходя из 
того, что искомое решение должно содержать какую-либо одну из дуг (1, /) или (7, 1), 
{, ЛЕК, не включенную в ЁЕ.. 
Если подмножество решений, определяемое точкой ветвления [, пусто, то оцен- 
ка ф, равна со. В этом случае хотя бы один из элементов 2; или Й, также равен со. 


Метод поиска оптимального решения задачи 


Алгоритм поиска у“ (К) строит оптимальное решение КЗСП на непустом мно- 


жестве её допустимых решений или устанавливает, что для данного графа Н не сущест- 
вует простых циклов, проходящих по всем рёбрам множества К, |К|>1. При |К| = 1 


КЗСП упрощается до задачи нахождения в графе Н кратчайшей цепи (р, ..., 4), 
соединяющей вершины {р, 4} Е К. Искомым решением задачи является цикл (р, ..., 4,р) 
стоимостью Гра + 4, где Гр — длина кратчайшей цепи (р, ..., 4). Известно, что её 


можно построить за время (т? ) [3]. 
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При изложении алгоритма поиска у“ (К) предполагается, что удаление строки 


(столбца) в исходной или сформированной матрице состоит в присвоении с каж- 
дому элементу этой строки (столбца). 


50. | 4, | — матрица стоимостей орграфа С = (Т, Е), |И =п, п>4, построен- 
Уп 


ного из связного графа Н = (Г, Ц), |Е| = 2} если {1, ЛЕ Ц, то (1, ]), (, )ЕЁЕ; граф Н не 
содержит мостов и точек сочленения, и степени всех его вершин больше 2; 4; е Я 


если (1, ЕЁ, иначе 4,=о, 2 — множество целых неотрицательных чисел; 


К -— непустое множество зафиксированных рёбер, образующих в Н паросочета- 
ние, |К| > 1; 

ИК) — множество вершин, инцидентных рёбрам К; 

О множество дуг, запрещающих образование циклов, которые не являются допус- 
тимыми решениями у(К), О=@; 

ИГ.) - расширение множества вершин И(К), ИГ.) = ©; 


положить О) = |4 | ; в дереве перебора матрице 2 соответствует корневая вершина; 
е п 


положить Аес = с; 

51. Для матрицы О выполнить модифицированный алгоритм Флойда-Уоршелла и 
получить в результате матрицу длин кратчайших путей М и матрицу маршрутизации И’. 

$2. Из матрицы М выделить рабочую подматрицу и преобразовать её в приве- 
денную рабочую матрицу Мьг, которая включает элементы матрицы М, принадлежа- 
щие строкам с индексами вершин из множества /„„С\ (ИСК) ИЕ.) и столбцами с ин- 
дексами вершин из множества [^^ ИК); [е» [с — множества индексов соответственно 
строк и столбцов матрицы М. 

$3. В приведенной рабочей матрице М» для всех элементов тра =0, р, ЧЕ 


= ИК) ИГ) найти оценку 
Й — 1 | 7 1 ' Й 1. ] М К цв М уй Й 
7 (р 4) Е Ш ое = ( ) ( ) (3) 


и определить максимальную из них 
Г(р, а) = тах (р, 4) | т =0] : 


по матрице маршрутизации И’ установить кратчайший путь (А, /) из вершины К 
в вершину /. 

54. Если {К, ПЕК, то множество решений разбивается на два подмножества: 
первое подмножество включает маршруты, проходящие по дуге (А, [), а второе — 
маршруты, проходящие по дуге (1, К); в дереве ветвления разбиение порождает две 
висячие вершины; 

получить из матрицы О две матрицы О’ и 0”: первую исключением строки Ки 
столбца [, а вторую — исключением строки / и столбца (; 

положить Д=р’и Р=р” соответственно для первого и второго подмножест- 
ва разбиения; 

применить для каждого подмножества разбиения правила формирования запре- 
щенных дуг О; 

в матрице О при каждом подмножестве положить равными оо все элементы, 
соответствующие запрещённым в нём дугам; 

выполнить шаги 51 и 52 для полученных матриц; 
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в подмножествах разбиения оценить стоимости решений по формуле (2); 

для обоих подмножеств выполнить проверку: если размерность матрицы Мах 
меньше 3 и оценка меньше Рес, то обновить значение Кес и запомнить соответ- 
ствующее решение как претендент на оптимальное, в противном случае, добавить 
висячую вершину к дереву ветвлений; 

перейти к шагу 57. 

$5. Если (А, ЕЁ, {К П®ВЮ, то множество решений включает подмножество 
маршрутов, содержащих дуту («, /), и подмножество, в которых её нет; 

определить для первого подмножества матрицу О’ следующим образом: если 
существует ребро {х, К} ЕК, а в матрице О при вершине ветвления элемент (х, А) не 
равен со, то получить путь (х, А, /[) и удалить из О строку х, столбец [, строку и стол- 
бец А; если существует ребро {1 у} ЕК и в матрице О элемент ([, у) не равен с, то 
получить путь (А, [, у) и удалить из О строку ^, столбец у, строку / и столбец [, если 
он не был удален ранее; удалить из О строку К и столбец [, если {х, К}, у} К; 

определить для второго подмножества матрицу О”, положив в матрице О 
элемент (А, /) равным оо; 

применить для первого подмножества правило формирования запрещенных дуг, в 
матрице О’ положить равным с элементы, соответствующие запрещенным дугам 0; 

положить О =р” и выполнить шаги 51 и 52; 

по формуле (2) определить нижнюю границу стоимости решений, входящих в 
первое подмножество; 

если размерность матрицы М» первого подмножества меньше 3 и его оценка 
меньше Кес, то обновить значение Кес и запомнить соответствующее решение как 
претендент на оптимальное, в противном случае, добавить висячую вершину к дере- 
ву перебора; 

применить для второго подмножества правило формирования запрещенных дуг; в 
матрице О” запретить значения элементов, соответствующих запрещенным дугам на со; 

положить О =” и выполнить шаги $51 и $2; 

по формуле (2) найти нижнюю оценку стоимости решений второго подмножества; 

если оценка меньше Кес, то присоединить висячую вершину к дереву перебора; 

перейти к шагу 57 ; 

56. (А, [) = (А, уььуо, ..., Ул У, ..., Унь Уь 0 = кратчайший путь из вершины / в 
вершину /; 

построить путь (а, (Е... ВИ где 


ра (х,К) ‚если {хе КиЧ =, в (1, у), если {1 у} Е Киа, =, 


©, иначе, ©, иначе; 


множество решений является разбиением на г + 2 подмножеств, в котором 
первое подмножество содержит маршруты, не включающие дуги (Км) пути 


(а, (Ель, ), В); каждый маршрут 5-го подмножества проходит по пути (©, 
[о )) и не проходит по дуге (у, 1,у,), в ("+2)-е подмножество входят марш- 
руты, включающие путь (анну, ), В | в дереве перебора добавляется г + 2 


висячих вершин; 
применить для каждого подмножества разбиения правило формирования 
запрещенных дуг О; 
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применить правило формирования множества (Г) и получить в результате для 
каждого 5-го подмножества расширение (Рф, и 5=2, "+1, где и (.) — рас- 
ширение в первом, и (и + 2)-м подмножестве разбиения; 

для каждого полученного подмножества найти матрицы ДО, 0..5, П[,, 
5=2, "+1, [= "+ 1; 1) чтобы получить матрицу Ду, элемент (Км) матрицы О при 
вершине ветвления полагается равным со; 2) матрица ),,, определяется удалением 


из О строк и столбцов, соответствующих промежуточным вершинам пути 


(а, (Ем... у В 3) матрица Д., 5=2, "+1, [=ги- 1 находится так: а) в мат- 


ай 
рице Д элемент (у,-1, у,) принимается равным со, 6) из р удаляются строки и столбцы, 
соответствующие промежуточным вершинам пути (а,(Е, му )); 


для всех полученных подмножеств выполнить проверку: если размерность мат- 
рицы Мь, меньше 3 и оценка меньше Кес, то обновить значение Кес и запомнить соот- 
ветствующее решение как претендент на оптимальное, иначе присоединить висячую 
вершину к дереву ветвлений; 

перейти к шагу 57; 

57. Если была обновлена граница Кес, то просмотреть все висячие вершины дере- 
ва ветвлений и удалить те из них, оценка которых больше или равна Кес; 

если дерево ветвлений не содержит висячих вершин, то процесс решения завер- 
шен, переход к шагу $38, иначе выбрать висячую вершину, имеющую наименьшую 
оценку и перейти к шагу $3; 

58. Если Кес равно со, то множество решений исходной задачи пусто, иначе опти- 
мальным решением задачи является решение со стоимостью Кес. 


Пример 


Для связного графа Н = ) и матрицы весов его рёбер [4; |. ‚ представ- 


* < 
ленных на рис. 2, требуется построить простой цикл у (К) минимальной стоимости, 
проходящий по рёбрам множества А={{1,2},{3,4}}, или установить, что множество 
маршрутов у(К) пусто. 


| 


©чмлмляияфчфь- 


5 || + © 8 | 
18 м1 8 © 
м8 м 62| 8 | а 315 
8 © 3% + 
юм юри 
8 © 8щшмы ам 
|8 мм 883% 
З < З © чм с 


Рисунок 2 — Граф Н КЗСП и матрица его стоимостей 


Граф Н не содержит точек сочленения, мостов и вершин, степени которых 
меньше 3. Имеем И(К)={1,2,3,4}, О =@ , У(Г.) =®. В корне дерева перебора Хо пола- 


гаем О = [4; |. ‚ Вее=о®. 
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Для матрицы Р модифицированный алгоритм Флойда-Уоршелла находит матри- 
цу кратчайших путей М и матрицу маршрутизации И’: 


И 3 О ВЕ. 1 о 3 4 5 6 7 8 
[|9 [6714] 5[3 [3 1 о [(.2) [а,53)[ а.) [а,5) ал.6[ ал) [а,3) 
оо [5:9 210.0 [ ® [(.3) [ 0.4) | (0,5) | ©.6) [0.6/7 0,3) 
6 [598 [25:4 |5 3(3,5,1)[ (8.2) [ ° [ (3,4) [ (,5) [| (3,6) [3,57 [ 6,8) 

М=4 |997 4 (4,1 | (42) | (43) [ ® | (4,5) (4,7,6)| (47) | (4,8) 
54529 94| 29 5 (5,1) | (5.2) | (5,3) | (5,4) | ® [(5,7,6) (5.7) | (5,8) 
6535 |4 ®|2 [10 6 (6,7.1)| (6.2) | (6,3) (6.7.4) (67,5) ® | (6,7) [(6,1,8) 
731514 9[2[2 [598 7| (2,1) [@.6.2)(,5,3)| @.4) | @,5) [| @,6) | ® |118) 
5 19210889 8| (8,1) | (8.2) | 8,3) [ 8,4) [ (8,5) 87,60 @Л) | ® 


Из матрицы М удалим строки и столбцы, соответствующие вершинам мно- 


жества Г\(И(Ю) ИЕ,)). В результате получим рабочую подматрицу [т | которая 


4? 
после приведения по строкам и столбцам преобразуется в приведенную рабочую 


матрицу Ми, = [т,. | р. 


[м„„ |. = 


8 © 1-14 
= 
Ц 


1 4 
оо 1 

9 0 |. 
6 0 

7 оо 


1-11 8 © 


3 
6 
7 
со 
3 


|-> 8 |= 
== З = 
фе о К ый К 


фр = 


По формуле (2) вычислим оценку Фх, . Сумма первых двух слагаемых равна 24, тре- 
тье слагаемое представлено суммой Д,> = пуп {и›, т} =0 и Дуд = пт { ид, т! =0, 
а четвёртое равно 0, поскольку Ё, =@. Следовательно Фх =25. 

В приведенной рабочей матрице Ме по формуле (3) находим для каждого элемента 
т’, =0 величину У(р,4).. В результате получим у (1,2)-1, у (1,3)=0, у (2/11, у (2,3)0, 
у (2.4)=0, у (3,4), у (4,3)1. Определим оценку К, /) = тах {у (ра) [т = 0} =1 и ©о- 
ответствующую ей дугу (1,2) орграфа С = (У,Е) ‚ {1,2} ЕК. 


Так как дуге (1,2) соответствует в графе Н ребро из К, то выполняются дейст- 
вия шага 54. Множество решений задачи разбивается на подмножество маршрутов, 
включающих дугу (1,2), и подмножество маршрутов, проходящих по дуге (2,1). В де- 
реве перебора разбиение означает появление двух висячих вершин Ху!, Х› и двух 


матриц 

13405067 8 213405 6 7 8 
ОО: 3: |997 1101717141613 |3 
3171090131215 |°%|5 31790931215 |°®|5 
473109 9|7 ь 417309 ®|9|7 
ан а 55 Ь [9 [9 [9 
66151090902 |<© 61315109020 < 
71310191212 °®|8 7 |1 91212 |8 
8131517908 81715171908] < 


Для матрицы О’ при вершине Л! Х\ О = {(2,1)}, поэтому элемент (2,1) принима- 
ет значение со; /(Ю)={1,2,3,4}, ИБЕФЯ. 

Модифицированный алгоритм Флойда-Уоршелла, входом которого является 
матрица О =р’, строит матрицы М и 
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К. 5-9. 25 26:96 98 1 3 4 5 6 Й 8 
ря 2 ® |063) бл 65 60 067 6 
3653 275 а [5 3 (3,50 [< | (3,4) [ (3,5) [ 6,6) [3,57 [ 6,8) 

мА 9 9 [|440 [а [| ® | (4,5) [4.7.6 | ал) [а 
54 [29 4[2[9} 95 (5,1) | (53) [64| ® [(5,1,6)[ 6/7) [ (5,8) 
655 и 4[%|2 | 10 6 (6.7.0 (6,3) [(6.1.4) (6/7,5)[ ® [| (6,7) [(6.18) 
[3141912 [2 [5 [8 71а [(4,5,3)[ 0.4) [ @,5) [ @,6) [| ® [4(5) 
83579108 ® 8| (8,1) [ (3) | 8/4) [ (8,5) @70[@Л | ® 


В матрице М при вершине Л! /„„={2,3,4,5,6,7,8}, [={1,3,4,5,6,7,8}. Рабочая 
подматрица матрицы М состоит из элементов строк с индексами вершин множества 
м (ИЮ®) ИЁГ)Е{2,3,4} и столбцов с индексами вершин множества /»/Г^\ ИК)={1,3,4}. 
После её приведения по строкам и столбцам получим матрицу 


Те В 24 не 
2101010 21901717 
Мия 3 ооо В = [6 [о [3 
41 [о |5 ЕЕ 


Найдём оценку фх,‚ вычисляемую по формуле (2). Сумма первых двух слага- 
емых в Фх равна 16, третье слагаемое 5. = тт{т.д,т/з} даёт 0, а четвёртое — вес 
4,2 =9 для дуги (1,2). Следовательно фх, =25. 


По матрице О” при вершине Х», для которой О={(1,2)}, И®)={1,2,3,4, ИБЕФХ, 
находим матрицы 


РЕ Е В. 2 3 4 5 6 7 8 
ео 6 [9 4.5.33 1 о [а,5.3) 4) | а.5) [а.7,6) ал) [а,8) 
3 [55] 3:12 [514 [5 3| (3,2) | < [ (3,4) | (3,5) [ 6.6) (6,57) 6,8) 

мА 9 [| 414.2 [43| в [| (4,5) [4.7.6] @7) [4,8) 
55 29% 429 ’” 5 (5,2) | (5,3) | (6.4) | ® [(5,7,6)| (5,7) | (5,8) | 
63514 2 [10 6 (6,2) | (6,3) (6.7.4) (6./7,5)| ° | (6,7) (6,1,8) 
75492] 2 [58 7(2,6,2) (1,5,3)| (4) | @,5) [@,6) [| с | (4,8) 
87151 7[9 |108 |® 8| (8.2) | (8,3) | (8.4) [ (8,5) [87,6 @7) | 


Поскольку /[„о„»С\ (ИК) И) { 1,3,4}, ис ИЮ®)={2,3,4}, то 


23-24 234 
11° |6 |7 110010 ]|1 
= =12+4 =1 НР: 
[т,, ]. 317 |°ю|3 › Мы 30 |© 01’ Фх, + +0+4> 6+9 5 
41713 |0 41010 [© 
Выберем вершину Х! точкой ветвления. В соответствующей ей матрице Ме 


оценки 7 (2,3)=0, у (2,4)=0, у (3,1)ЕТ, у(3,4)=0, у (4,3)=1, Г(З,П)=1. Ветвление иници- 
ирует элемент (3,1), который в матрице маршрутизации И’ представляет путь (3,5,1). 
Поэтому выполняются действия шага 56. 

Так как {3,4} Е А, 443 © ‚то а =(4,3), но В =© ‚, поскольку ребро {1,2} е К имеет в 
матрице ШО” вес 412=®. Получен путь (4,3,5,1), который вызывает разбиение множест- 
ва решений, представленного вершиной АХ! на три подмножества Хз, Ха, Х5. Каждый 
маршрут в Хз не проходит по дуге (3,5), в Х4 содержатся все маршруты, включающие 
дуги (4,3),(3,5), а в Х5 включены маршруты, проходящие по дугам (4,3), (3,5),(5,1). В де- 
рево перебора добавляются три вершины, исходящие из точки ветвления Х\ (рис. 3). 

Для множества маршрутов Хз, не включающих дугу (3,5) О={(2,1)}, ИК) = 
— {1,2,3,43, ИБ @ , @35=00 , @1=00 , 
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В: ВЕ Е ВЕ 1`.3. 245-10. 9.8 
2109017171513 |°|[7 2191717151315 17 
3171%°|3 19015 |9 | 5 3171031915171 5 
о ЕЕ ИЕН 
51412191%9191219Р 5141219104129 [? 
66509 ®|2 |< 615151111 41°%| 2110 
7131090191212 |©®©|[8 713149212 ©®1[8 
8131579585 831517] 9 [10] 8 © 
1 3 4 5 6 у 8 
21 ох | (2,3) | (2,4) || (2,5) | (2,6) |(2,6/7)| (2,8) 
3| (3,1) | <® | (3,4) |(3,6,7,5) | (3,6) |(3,6,7)| (3,8) 134 
р (4,1) | (4,3) | © (4,5) (47,6) | (47) | (458) у 2 00| р 5. 
5| (5,1) | (5,3) | (5,4) ®  |(5,7,6)| (5.7) | (5,8) "30 |ю|0|’`*” 
6 |(6,7,1)| (6,3) |(6,7,4)| (6,7,5) | ® | (6,7) |(6,7,8) 41010 |< 
7| (7,1) |(1,5,3)| (7,4) | (7,5) | (1,6) | ® | (1,8) 
8 (8,1) | (8,3) | (8,4) | (8,5) |(8,7,6)| (87) | 


Множество маршрутов Ха, проходящих по дугам (4,3),(3,5) и не включающих 
дугу (5,1), формирует О={(2,1),(5,4)}, ИЮ)={1,2,3,4}, ИД-{5}, 2 00, 454 00, Ч 00, 


14678 14678 ве 4.28 
257 [3 [7] 2 7357] 2[ о [4 [©.6 ©67 © Е. 
Сео [9 [29| 515 [94 [29| _ 5[5.7.0[ (5.4) 5.7.6) (6.7) [ (5.8) 
в ооо [2 о М5 Ио пб! 56979 5 [6 "= таг. 
392 м8 ТЗ 9208 700104 |@0| о |(@9 
Я © 8 © 83171108 © 8 (81) (84) (876) ®7) | © 


После приведения матрицы [т,, |, получим матрицу 


2 

1 
Ми=2 [© 
510 


8©- 


Для множества маршрутов Х5, проходящих по дугам (4,3),(3,5),(5,1), О={(2,1)}$, 
июк)= { 1 ,2,3,4} ‚ Ир © > 424=00 ‚ 


4678 Е 4 6 й. 8 
27 3®|7 27357 210.4) | 0.6) 0.67) | 0,8) 4 
рб [оо [2 [5 |, 6 5 [2 [10], = 6 (6.7.4) ® | (6,7) [(6,7,8)|, Мы= 2[0] 
92 8179258 7! (1.4) | (0.6) | ® [| (18) 
87 [| 8 [© 8171 10| 8 [с 81 (8,4) [(8.7.6)| 87) [ ® 


Рабочая подматрица Мю состоит из одного элемента (2,4), отсюда следует, 
что получено допустимое решение задачи. В матрице И’ элемент (2,4) представлен 
дугой (2,4) стоимостью фл= 7. Стоимость решения, содержащего дуги (1,2), (2,4), 
(4,3), (3,5), (5,1), равна нижней границе Фх; = 412 аз | 435 51 | фа = 9+3+2-+ 


+47 = 25. Поскольку значение Кес > фх.,то Кес =фх. =25. Оценка фу, является 


наименьшей из всех оценок висячих вершин дерева перебора. Следовательно у*(А)= 
= (1,2,4,3,5,1), СУЖ(К))25. 
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Модификация метода Литтла для решения кольцевой задачи... 2М 


Рисунок 3 — Дерево перебора, построенное для данных примера 


Выводы 


В статье сформулирована ограниченная версия известной задачи о сельском 


почтальоне — кольцевая задача о сельском почтальоне. Предложен метод её решения, 
развивающий классический алгоритм ветвей и границ для решения задачи ком- 
мивояжера (алгоритм Литтла). Разработанный метод также может быть применён 
для решения гамильтоновой задачи о сельском почтальоне и гамильтоновой задачи 
коммивояжера. 
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Модифкац методу Литла для розв’язання к«льцевот задач! про сльського листоношу 

У статт! приведено формулювання 1 анал!з узагальнення гампльтоново! задач! про слльського листоношу. 
Пропонуеться точний метод Й розв’язання, який розвивае класичний алгоритм Литла. 
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